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We present a simplified version of the atomic dark matter scenario, in which charged dark con- 
stituents are bound into atoms analogous to hydrogen by a massless hidden sector U(l) gauge 
interaction. Previous studies have assumed that interactions between the dark sector and the stan- 
dard model are mediated by a second, massive Z' gauge boson, but here we consider the case where 
only a massless 7' kinetically mixes with the standard model hypercharge and thereby mediates di- 
rect detection. This is therefore the simplest atomic dark matter model that has direct interactions 
with the standard model, arising from the small electric charge for the dark constituents induced by 
the kinetic mixing. We map out the parameter space that is consistent with cosmological constraints 
and direct searches, assuming that some unspecified mechanism creates the asymmetry that gives 
the right abundance, since the dark matter cannot be a thermal relic in this scenario. In the special 
case where the dark "electron" and "proton" are degenerate in mass, inelastic hyperfine transitions 
can explain the CoGeNT excess events. In the more general case, elastic transitions dominate, and 
can be close to current direct detection limits over a wide range of masses. 
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In recent years there has been increased interest in dark 
matter models in which the dark sector has some of the 
richness of the visible sector, such as hidden gauge inter- 
actions [3-0 and flavor. A natural possibility to consider 
is an unbroken U(l) gauge symmetry that would give rise 
to bound states, i.e., atomic dark matter |8|-|ll|. In this 
case there should be a further resemblance to the visible 
sector in that the dark matter must be asymmetric [l2j 
in order to have the right abundance; otherwise the U(l) 
coupling must be so weak that recombination in the dark 
sector does not occur efficiently and the would-be atoms 
remain predominantly ionized fl3j . 

Atomic dark matter can have interesting properties 
with respect to direct detection, because of the possibility 
of inelastic scattering to excited states of the atom, no- 
tably through hyperfine transitions. In previous studies 
it has been shown that inelastic transitions can help to 
reconcile the DM interpretation of CoGeNT even ts | 14j 
with null results from XenonlO [l5j]. In refs. (T^, it 
was assumed that the hyperfine transitions were medi- 
ated by the kinetic mixing of the photon with a massive 
dark vector boson that couples to the axial vector cur- 
rent of the DM. In the present work we explore a simpler 
possibility @: one can rely upon mixing of the photon 
with the massless 7' that is already present due to the un- 
broken U(l) gauge symmetry. This is a very economical 
model, while still endowed with the rich phenomenology 
of the atomic DM scenario. 

Because of the kinetic mixing, the constituents of the 
dark atoms acquire small electric charges ee; thus the 
interactions that give rise to direct detection are electro- 
magnetic. We show that direct searches in fact give the 
strongest bounds on e in this model; thus such detections 
could be imminent. In fact in the special case where the 
two constituents have equal mass, we offer an interpreta- 
tion for the excess events reported by CoGeNT, relying 
upon inelastic interactions. 



We follow the notation of ref. [10( by denoting the dark 
analogues of the proton, electron, and hydrogen atom by 
p, e, H. The Lagrangian is 
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where F is the electromagnetic field strength and F' is 
that of the massless 7'. D' is the covariant derivative 
with respect to F': D' = d ± igA', where g is the U(l)<z 
dark coupling constant. We ignore the small mixing of 
the 7' with the Z boson, and hence refer to F as the elec- 
tromagnetic rather than the hypercharge field strength, 
e is the gauge kinetic mixing parameter. The gauge bo- 
son kinetic terms can be diagonalized to first order in 
e by letting F' = F' + IF. Then A'^ couples only to 
the dark current q_J% while the photon A^ couples to 
eJ em + cg J d E G3- The DM particles thus acquire 
millicharges eg = ee under the electromagnetic U(l). We 
will refer to e rather than e in the remainder of the paper. 



2.1. m e >C m p case 

Dark atom interactions. The low-energy interac- 
tions of H are screened due to its net charge neutrality. 
Let us first consider the generic regime where m e <C m p . 
In this limit, the Fourier transform of the H electric 
charge density is given by /?h — ee a' 2 q 2 /2 at low wave- 
number q <C l/flf), where a' = l/(o/m e ) is the Bohr 
radius of H, and a' = g 2 /4tt. We assume that m e -C m p 
in our approximation for a' . The factor of q 2 in p cancels 
the factor of 1/q 2 coming from the gauge boson propaga- 
tor in Coulomb gauge, so that the scattering of H on a 
proton in a DM detector will not be long-range, but will 
instead appear as a contact interaction, with 
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where /i is the reduced mass of the p-H system. This 
expression relies upon the Born approximation, whose 
validity depends upon the properties of the central po- 
tential experienced by p due to H: 

The Born approximation is justified if |V(a[))| <C 
l/(/xa' ), which implies 

e — <^ 7 Ih. - m e K + TO H) / 4 x 

a' /i TOp ran 

We will see that this condition is satisfied for param- 
eters of interest for direct detection. In comparing a p 
to direct detection bounds, we must take into account 
that H interacts only with protons and not all nucleons. 
This weakens the experimental limit on a p by a factor of 
(Z/A) 2 for a target with atomic number and weight Z, A. 
For Xenon, {Z/A) 2 = 0.17. We define a p . oS = (Z/A) 2 a p 
to facilitate comparison with the Xenon excluded region. 

The cross section for direct detection is proportional 
to the combination 

e 2 (l + ze) 4 _ 

when we reexpress fi and a' in terms of m e and mm, 
where x e = m e /m p = m e /mn. The values of ft that 
are interesting for direct detection can be read from fig. 
[TJ where contours of constant /3 (labeled by the value 
of log 10 /3), are plotted in the mn-a Pie g plane, on top of 
constraints from the XenonlOO experiment. To give more 
concrete examples, let us take a' = 0.1 and x e = 0.1, 
which tend to give a small ionization fraction and there- 
fore are more robust with respect to structure formation 
constraints. (The ionization fraction / is computed in 
analogy to that of visible hydrogen in ref . , and goes 
roughly as / ~ 10~ 10 a'~ 4 m e m p GeV~ 2 .) The values of e 
needed to saturate the XenonlOO bound [l8[ for a range 
of ran are shown in figure [2j These scale as x 2 a' 2 for 
different values of x e and a' . The values shown in fig. [2] 
easily satisfy condition 

In principle, inelastic scattering can also occur, in 
which the internal atomic state changes. The lowest en- 
ergy excitation available is the hyperfine transition, with 
energy gap AE = ^a' 4 m 2 /m p . For our fiducial param- 
eters a' — x e =0.1, this would correspond to ~ 30 keV 
if mil — 10 GeV. Although this may be an interesting 
value for direct detection, the rate of such transitions is 
suppressed compared to the elastic ones by a factor of 
a' 4 x 2 (mn/m p ) 2 so they make a subdominant contribu- 
tion if mn ^ TeV. Transitions that excite the e orbital 
state are less suppressed by powers of a', but have a 
larger energy gap, and so are also irrelevant. 

Dark ion interactions. An interesting feature of 
atomic dark matter is that a small fraction remains in 
the ionized state, which has a larger cross section on nu- 
cleons than does the atomic state, because its charge is 
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Figure 1: Diagonal lines: contours of constant B, eq. (|5]l. in 
the m,H-o"p, e ff plane. Lines are labeled by the value of logi n 8. 
Background shows limits from the XenonlOO experiment [l8| ] . 
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Figure 2: Lowest curve: values of e that saturate the 
XenonlOO direct detection bound, as a function of run, for 
elastic scattering of atoms in the model with a! = m e /m D = 
0.1. Upper curves: upper limits on e from CMB (ref. (lap, 
capability of p ions to penetrate 1 km of rock (this work), 
CoGeNT limit (adapted from [l^]) and XenonlOO limit (our 
estimate) from detection of p ions. Note that the ion detec- 
tion limits may not apply (see section 2.1). 

unscreened. Therefore one might question whether the 
scattering of ionic DM could dominate the direct detec- 
tion signal, even though it is a subdominant component. 
Ref. 19] has argued that the ionized component would 
necessarily have been blown out of the galactic disk by 
supernova shock waves 20] if e was in the range required 
for direct detection. The galactic magnetic field subse- 
quently shields the disk from being repopulated by the 
millicharged ions, unless [2l| the strong scatterings be- 
tween the ions themselves sufficiently randomize their 
directions contrary to the Lorentz force from the mag- 
netic field. Adapting the estimate (8) of [2l| for the dis- 
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tance scale I, over which randomization occurs, to our 
models with a 1 — 0.1, m e /m p = 0.1 and using the wh- 
dependent ionization fraction from fig. 1 of ref. [1(3], we 
find that I = 300 pc, which is larger than the height 
of the galactic disk and thus ineffective for overcoming 
magnetic shielding. 

Nevertheless, in case there may be some other way of 
evading this argument, we indicate by the dashed line of 
fig. [2] the value of e at which dark ionic scattering would 
saturate the CoGeNT signal, which is adapted from ref. 
[l9| by taking into account the ionization fraction / men- 
tioned above. Thus even if the ions penetrate to the 
earth, their signal is much weaker than that of the atoms 
unless ran > 100 GeV. For consistency one must also 
check that dark ions can penetrate ~ 1 km of rock, which 
we have done along the lines of refs. 

GLU- Fig. ® shows 
that much larger values of e are needed to stop the ions 
in the earth than to detect them. 



2.2. Special case m e — m p 

If for some reason (e.g., a discrete symmetry) m e = 
m p , the matrix element for elastic scattering vanishes in 
the Born approximation since the average charge density 
in the H atom vanishes. This is an interesting situation 
since then inelastic scattering can be the dominant effect 
for direct detection. The hyperfmc splitting is given by 
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assuming gyromagnetic ratios g e = g p = 2 and mu = 
2m p . The transitions from the spin singlet to triplet 
atomic DM states are dominated by the spin-orbit cou- 
pling between the proton and the constituents of H, 
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where fl e — ga e /m e is the dark magnetic moment of e 
(hence efl e = eea e /m e is its normal magnetic moment) 
and r is the distance between p and e. (Notice there is no 
reduction by 1/2 for Thomas precession since the electron 
rest frame is effectively inertial.) We neglect the spin- 
spin couplings because these give rise to spin-dependent 
interactions with the nucleus that are suppressed due to 
the lack of coherence. 

The squared matrix element, summed over the final 
spin states of the triplet, and taking into account the 
equal contributions from e and p, is given by 
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where p, p' are the initial and final momenta of the pro- 
ton (in the rest frame of H), q = p — p' is the momentum 



transfer, C = ee 2 /m e , and (i p /p = m p mji/[(r 



mu)p] 



is from normalizing the incoming plane wave to unit cur- 
rent density [22| ■ From this we obtain the inelastic differ- 
ential cross section for protons on dark atoms, which can 
be rescaled to represent the nucleus-atom cross section 
by including the factor of Z 2 to sum over the individ- 
ual proton contributions, letting fi p — > /ijv, and inserting 
the Helm form factor F^(q) to account for the nuclear 
structure: 



ln~ m\ Vp ] ql 
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where now p, p' stand for the initial and final momenta 
of H (in the lab frame), q = p — p' as before, and we used 
rn e = |m H . 

Let us compare to the corresponding result (40, 42) of 
[Io| in the m e = m p limit: 



dcr N 4a a /j, N q 2 
dQ ~ tt m 2 p / e 4 ff H 
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By equating © with (fT0|) using the preferred values for 
fitting to the CoGeNT data, A = 73, q = y/2m N E R - 26 
MeV for recoil energy Er — 5 keV, ti~2x 10 _3 c, and 
f~ s 2 - 10- 30 cm 2 [TjJ], we find that 
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to explain CoGeNT. The splitting E hi = 15 keV and 
atomic mass mn = 6 GeV imply a' = 0.062. We confirm 
the estimate (fTTj) by comparing a p ~ 7r[8e a n P v / (q toh)] 2 
(c.f. eq. ©) to the determination a p = 10~ 38 - 3 cm 2 from 
fig. 12 of ref. [23|], which was the first to propose inelastic 
scattering (also with a mass difference of 15 keV) as an 
explanation for the CoGeNT observations. We note that 
dark ions will not penetrate 1 km of rock for such large 
e, considering fig. [2j 

3. OTHER CONSTRAINTS 

Laboratory and supernova bounds. Unlike mod- 
els in which the 7' has a mass, in ours no coupling of the 
7' to visible sector matter is induced by the kinetic mix- 
ing. Therefore a variety of bounds that would pertain 
to massive 7's, from beam-dump experiments, contribu- 
tions to the anomalous magnetic dipole moments of the 
electron and muon, and supernova emission of 7', do not 
apply here. Note that the dark matter is too heavy to 
be in equilibrium in supernovae. Accelerator constraints 
for millicharged particles have a large open window for 
e < 0.1 and masses > 1 GeV of interest here (24l |. 

Exotic isotopes. Millicharged dark ions with ioniza- 
tion fraction / = 10~ 4 would ostensibly have a flux of 
2 x 10 20 /s on the earth, and they would bind to nor- 
mal nuclei unless they are unstable against thermal fluc- 
tuations, requiring e < 10" 3 [1, H)]. With e ~ 10~ 2 , 
they would be stopped in ~ 1 m of the atmosphere (10 44 
atoms) and produce a relative abundance 10~ 7 of exotic 
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isotopes over 10 Gyr. In the mass range covering m e = 3 
GeV, heavy isotope searches have excluded abundances 
of 10~ 18 - 5 for deuterium from D 2 H and 10~ 14 for he- 
lium [13. 

However, there are a number of reasons why these lim- 
its would not apply to our model. We evade the first 
one because e binds much more strongly (400 eV) to 
oxygen than to deuterium (3 eV) for e = 0.01, making 
D[De]0 highly unstable to decay into D2[Oe]. The sec- 
ond is ameliorated by realizing that He has a lifetime of 
t = 10 6 y in the atmosphere [28], reducing the estimated 
abundance by a factor of r/(10 Gyr) to 10 -11 . This does 
not yet take into account the magnetosphere which very 
effectively shields the earth from slow charged particles, 
including 3 GeV dark ions with e ~ 10~ 2 , whose gyrora- 
dius at the top of the atmosphere is ~ 0.01 earth radii. 
Solar x-rays are sufficiently energetic to break up the He- 
e bound state (binding energy 5 eV) and allow e to rebind 
much more strongly to N or O in the atmosphere. Even 
though ref. [20] conservatively limits e < 0.005 for super- 
novae to be able to efficiently expel 3 GeV ions from the 
galaxy, which is marginally smaller than our preferred 
value, the galactic B field does prevent new ions from 
entering the galaxy, and this could decrease the expected 
flux. Moreover our determination of e could decrease by 
a factor of 10 5 in light of CoGeNT's recent reanaly- 
sis of their background events [29(. Thus it is far from 
clear that exotic isotope searches rule out our e ~ 10~ 2 , 
m e = TO p = 3 GeV model. 

Big bang nucleosynthesis (BBN). The current 
bound of no more than one additional neutrino [30| in 
the plasma at BBN gives weak constraints on e since the 
extra 7' counts as only 8/7 of a neutrino. Thus even if 
7' remains completely in equilibrium, it only exceeds the 
95% c.l. bound by 0.14 of an additional neutrino species. 
As expected, the constraints on e are quite weak from 
demanding decoupling of processes keeping 7' in equi- 
librium early enough to dilute this small fraction of a 
species. 

Neutron stars. It has been pointed out that asym- 
metric bosonic dark matter must interact extremely 
weakly with baryons in order to avoid the destruction 
of neutron stars due to accumulation leading to gravita- 
tional collapse [3ll433| . This restricts the cross section to 
values below that needed for direct detection, but it re- 
lies upon the absence of Fermi pressure for bosonic DM. 
Once the atoms become so dense that they are over- 
lapping, their constituents start to behave as a Fermi 
gas, and thus avoid collapse until much higher densi- 
ties, as long as the overlap of atoms occurs before they 
are within their Schwarzschild radius, R s . For N dark 
atoms of mass TOh, R s — 2Nmn/Mp and the average 
separation between atoms is Ar = Rs/N 1 / 3 . The crite- 
rion to avoid collapse is that Ar <C &o at the bosonic 
Chandrasekhar limit N ~ (M p /ra H ) 2 [31]. We find 
that Ar/aQ ~ a'm e / (m^ 3 Mp^ 3 ) -C 1, so the system 
is no longer atomic, but a plasma of dark ions. In 
this case the much weaker fermionic Chandrasekhar limit 



N ~ Mp/(me'' 4 mp 4 ) [34| applies, and there is no con- 
straint from neutron stars. (Ar/a' Q ~ a' ^Jm e /m p is still 
sclf-consistently < 1 in this case.) 

Halo shape constraints. The most serious con- 
straints on atomic dark matter self-interactions come 
from their distortions of the shapes of elliptical DM cores 
of galactic clusters. To avoid relaxation to more spher- 
ical shapes, ref. [35| finds the constraint a /ran < 0.02 
cm 2 /g on the cross section for H-H collisions. Ref. [Toj 
argues that a = 47t(k<2q) 2 with 3 < n < 10. Together 
with the bound on a, this would rule out our inelastic 
model with toh = 6 GeV. However, assuming that the 
elastic cross section for H-H scattering computed in [36[ 
is the same as for H-H, we find that k = 0.16 at the 
relevant energy E ~ (v/a 1 ) 2 ~ 10~ 3 in atomic units. 
The ellipticity constraint is then satisfied for the inelas- 
tic model, with wh ^ 2 GeV. For the elastic model with 
a' = m e /m p = 0.1, it implies mn ^ 4 GeV. Bounds on 
a from the Bullet Cluster give weaker constraints [37l ]. 

Cosmic microwave background (CMB). Bounds 
on e were obtained on millicharged DM in ref. [l|| con- 
sidering a variety of physical processes; the most strin- 
gent limits were obtained from demanding that dark mat- 
ter has decoupled from the photon-baryon plasma before 
recombination, under the assumption that the DM was 
fully ionized. These are weakened in our model due to the 
screened electromagnetic interactions of the atoms or the 
ionization fraction / being small. By extrapolating the 
results of fig. 1 of ref. [l(| to m e = m p = 3 GeV we find 
that / = 10 -4 at a' = 0.062. This is well below the limit 
/ = fiion/^atom < 0.007/0.11 from distortions of the cos- 
mic microwave background [38]. One also requires that 
the DM be out of kinetic equilibrium with the baryon- 
photon plasma before recombination. For the atoms, the 
dominant interaction is 7H — > 7H through the Compton 
cross section a c = 327re 2 a 2 /3m 2 . We find the weak limit 
e < 0.02, which is marginally consistent with (TTT|) . Thus 
we evade the strongest constraints on e in fig. 1 of [19j in 
our inelastic scattering model with m e = m p . 

In the case m e < m p where elastic scattering domi- 
nates in direct detection, we need not consider how much 
the bounds of ref. [l9| are softened by having dark atoms 
rather than ions. In this case, the XenonlOO limit on e, 
at least for the model a' = rrie/m p — 0.1, is well below 
the most stringent limit of I19j ] even if the dark matter 
was fully ionized. This is shown in fig. [2j 

Virialization of dark matter. The next-strongest 
bounds in [l!| arise from heating of the DM by scatter- 
ing from baryons, which would interfere with DM col- 
lapse and virialization in galaxy formation. However in 
the m e = rn p case, the tree- level exchange of photons 
between H and normal H atoms vanishes because of the 
perfect mutual screening of e and p. Thus there is no ef- 
ficient way of heating up the dark atoms; only the small 
ionized fraction suffers this fate, and it should have a 
negligible effect on galaxy formation. For the m e <C m p 
case, our XenonlOO bound on e is more restrictive than 
the virialization bound. 
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4. CONCLUSIONS 

Wc have reconsidered a minimal alternative for atomic 
dark matter, namely that at the Lagrangian level its ionic 
constituents couple only to one massless 7' gauge boson, 
responsible for their binding, but due to kinetic mixing 
of the 7' with the photon, they acquire small charges 
±ee giving the dark atoms a weak coupling to the ordi- 
nary photon. We have shown that e is sufficiently un- 
constrained so that the first evidence of this interaction 
might come through direct detection of the dark atom. 

There are two interesting regions of parameter space 
for the model with respect to direct detection. If the 
dark "electron" is much lighter than the dark "proton," 
m e -C m p , then dark atoms scatter primarily elastically 
on nuclei, and the cross section can be close to the limits 
from direct detection over a wide range of masses run- 
However if m e = m p , elastic scattering is suppressed and 
inelastic hyperfine transitions dominate. We find that if 
the 7' gauge coupling is a' = 0.06 and m e = m p = 3 
GeV, the hyperfine splitting is ~ 15 keV, and e = 10~ 2 
gives the right cross section for explaining candidate DM 
events reported by CoGeNT. 



This model has further testable implications. The de- 
tailed effects of atomic DM on the CMB, depending upon 
exactly how and when it recombines, have yet to be stud- 
ied. An important unexplored aspect is the fraction of 
atoms that combine into molecular H2. This could re- 
sult in multiple signals for direct detection, where dark 
atoms or molecules interact with nuclei, producing re- 
coils at energies corresponding to several different DM 
masses. Searches for exotic isotopes consisting of bound 
states of e and normal atoms appear to be tantalizingly 
close to ruling out the e = 0.01 model, but more work 
must be done to quantify the predicted abundances. The 
framework appears to be consistent with strong DM halo 
ellipticity constraints, but an independent determination 
of the H-H scattering cross section at the relevant en- 
ergies should be done to confirm this. The origin of the 
DM asymmetry [Tl| should also be addressed. We hope 
to return to these issues in the near future. 
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